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Abstract. We obtain high-resolution echelle spectra of evolved Hipparcos 
stars within 100 parsec for detailed abundance analysis. We are particularly 
interested in subgiants because their ages can be derived accurately. Here, we 
present preliminary results for 10 subgiants and discuss the future work for this 
project. Our main scientific interest is to incorporate our results to constrain 
the star formation history in the solar neighborhood. 



1. Introduction 

Accurate parallax measurements of solar neighborhood stars from the Hippar- 
cos catalog have enabled detailed studies of stellar properties including stud- 
ies of the star forma tion history (|Hernandez. Valls-Gabaud. Gilmord l2000l : 
iBertelli h Nasil Eboih By looking at stellar abundances at different ages, we 
can constrain the star formation history of the solar neighborhood and study 
the history of chemical evolution. Subgiants are of particular interest because 
they lie in an area where isochrones are well separated, enabling more accurate 
age determination. 

We have initiated a project to obtain high-resolution echelle spectra of 
evolved Hipparcos stars in a volume-limited sample within 100 parsec. We are 
performing detailed abundance analysis to investigate the chemical evolution 
and the star formation history of the solar neighborhood. 

Currently, we have obtained spectra for about 50% of our sample. We 
are attempting automated analysis to obtain fundamental attributes such as 
[Fe/H], effective temperature, surface gravity, and age. Here, we present initial 
results for 10 subgiants over a range of effective temperatures. We compare our 
results using several different methods to literature values, and discuss future 
development for this project. 
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2. Sample 




- v 



Figure 1. (a) Color-magnitude diagram showing all Hipparcos stars within 
100 parsec (small points) with sample cut lines (solid lines). Ten subgiants 
presented here are shown as large points, (b) Number of observed stars to 
date as a function of color. 



We select evolved Hipparcos stars with acceptable parallax uncertainties 
in the Northern Hemisphere within 100 parsec that have My < 4 and B-V > 
0.38, as shown in Fig. ^a). To probe histories of star formation and chemical 
evolution, abundance analysis of subgiants is ideal because we can use the well- 
separated isochrones of the subgiant branch to obtain accurate ages. Therefore 
we start our analysis on subgiants first, of which the 10 that we are presenting 
here are also shown in Fig. ^a). They are chosen to cover a wide range of 
effective temperatures to test the accuracy of our automated method. 

We are randomly selecting about 500 stars to observe as a representative 
sample. We exclude known binaries and multiple-star systems to avoid com- 
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plications such as abundance ambiguity in spectroscopic binaries. We will also 
exclude stars with high rotational velocities and significant on-going mass loss. 

By collecting our own data, we ensure a uniform sample, consistency in 
data reduction, and sufficient SNR to obtain abundances of various individual 
elements. We obtain high-resolution echelle spectra with the ARC 3.5-m tele- 
scope at the Apache Point Observatory, and use standard techniques for spec- 
trum extraction and calibration. To date, we have obtained over 300 spectra, 
tallied in Fig. ^b), and are performing initial abundance analysis. 

3. Analysis 



Fe I 




5813 5814 5815 5816 



Figure 2. An example of automatic nonlinear gaussian fitting for an Fe i 
line. The best fit is shown as solid line with flat continuum. The other solid 
line is the spectrum of HIP43454. Value of W shown is the equivalent width 
of the fit. 

Due to the number of stars in our sample, we are attempting to automate 
the analysis wherever possible. We measure equivalent widths (W) of selected 
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"clean" Fe i and Fe n lines using nonlinear least squares fitting with a gaussian 
profile. An example is shown in Fig. [21 For the 10 stars discussed here, we also 
manually measure W for comparison with t he automatic results. 

We use these measurements in MOOG (Hneden|[2002) with a grid o f lKurucd 
(1993) model atmospheres based on initial estimates of stellar properties. Our 
grid is interpolated from existing model atmospheres and has effective temper- 
ature (T e ff), logarithmic surface gravity (log g), and mircoturbulence velocity 
(£) steps of 50 K, 0.1 cm s~ 2 , and 0.1 km s _1 respectively. For each star, we 
obtain [Fe/H], T q q, log g, and £ by finding zero slopes of Fe i abundance versus 
excitation potential (EP) and reduced equivalent width (RW), and ionization 
balance between Fe i and Fe n. We can compare the log g obtained from parallax 
and T g g from color with those from spectroscopy. 

With a derived metallicity, we fit Pad ova isochrones ( Girardi et al.l f2002) 
through a star to obtain its mass (M), age (r), and "trigonometric" log g. These 
isochrones are interpolated for sufficient fitting via closest-point match in the 
L,T e g-plane and have [Fe/H] intervals of 0.05 dex. For the stell ar luminosities, 
we ap ply corrections to the absolute magnitudes as suggested bv lLutz fc Kelkerl 
(1923). '. n future work, we will improve isochrone fitting by including probability 
weighting to allow for uncertainties in the photometry and distances, and overlap 
between isochrones. 



4. Initial Results - 10 Subgiants 

We present preliminary results for 10 subgiants. They are selected from the 
color-magnitude diagram to span a range of temperatures in order to evaluate 
the accuracy of our method over varying stellar properties (in this case, T e ff)- 
We perform equivalent width measurement, abundance analysis, and isochrone 
fitting as described in the Analysis section. In Tabled we tabulate the results 
from the automated W fitting. 



Table 1. Results for the 10 subgiants via automatic W measurements. 



HIP 


^eff M 


log g [cgs] 


e [km/s] 


[Fe/H] 


M [Mo] 


r [Gyr] 




± 50 


± 0.1 


± 0.1 








13679 


6100 


3.5 


2.0 


-0.08 ± 0.05 


1.35 ± 0.22 


3.2 ± 2.0 


6869 


6000 


3.5 


2.0 


-0.02 ± 0.09 


1.22 ± 0.12 


4.2 ± 1.5 


9726 


5950 


3.9 


1.8 


-0.15 ± 0.07 


1.17 ± 0.08 


5.1 ± 1.8 


7276 


5850 


4.1 


1.6 


+0.05 ± 0.05 


1.21 ± 0.03 


4.9 ± 0.4 


43454 


5750 


3.8 


1.4 


-0.01 ± 0.05 


1.23 ± 0.06 


4.8 ± 0.7 


8159 


5650 


4.1 


1.3 


+0.05 ± 0.04 


1.12 ± 0.03 


7.1 ± 0.5 


39681 


5250 


3.5 


0.9 


-0.43 ± 0.05 


1.35 ± 0.69 


2.8 ± > 5 


52316 


5150 


3.9 


1.0 


+0.01 ± 0.06 


1.27 ± 0.05 


4.4 ± 0.5 


58080 


5100 


3.7 


1.1 


+0.14 ± 0.07 


1.55 ± 0.25 


2.5 ± 1.7 


61995 


4800 


3.5 


0.9 


+0.26 ± 0.11 


1.12 ± 0.39 


8.0 ± > 5 



HIP13679 lies at the hook of main sequence turn-off, HIP39681 has consid- 
erable parallax uncertainty (— = 0.16), and HIP61995 turns out to be on the 
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red giant branch. These cause substantial uncertainties in Ai and r values. In 
the future, we will refine our selection criteria to minimize these issues. 




Figure 3. Deviation of our [Fe/H] results from literature values for (a) au- 
tomatic and (b) manual W measurements. Error bars for [Fe/H] include 
both our uncertainties and assumed 0.1 dex errors for literature values. Error 
bars for T e g is from model atmosphere grid spacing of 50 K. HIP58080 and 
H1P61995 do not have published values for comparison. 



We compare our iron a bundance from automated and manual W measure- 

ment s with literature values ([Cavrel de Strobel et al.ll997l : ICavrel de Strobel. Soubiran. fc E alitel 
120011 : iNordstrom et al.l 12004) in Fig. 131 On average, the automated measure- 
ments agree better with published values, except for HIP 13679. More data on 
the lower and higher T q q ends are required to determine if there is any [Fe/H] 
trend with effective temperature. 

We compare the derived stellar parameters from th e automatic and manual 
W m easurements in Fig. @J We use T eff from V — I dBessel. Castelli. fc Pled 
1998) for visual convenience. The "automatic" method underestimates [Fe/H] 
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Figure 4. Deviation of manual W measurement results from the automatic 
for (a) metallicity, (b) effective temperature, (c) surface gravity, and (d) mi- 
croturbulent velocity. Effective temperatures derived from V — I for the 10 
subgiants are plotted on the X-axis. 



and overestimates £ compared to the "manual." However the results for log g 
and T q q are comparable except at > 6000 A. 

Two possible reasons for these deviations are systematic differences in the 
two W measurement methods and slopes affected by outliers at the upper or 
lower limits of EP and RW during abundance analysis. This calls for more data 
and careful visual inspection of the results. For the former, we could fine-tune 
the automation, such as implementing dynamic continuum determination and 
line templates as a function of temperature. For the latter, we could add more 
lines near the EP and RW edges, or formulate an analysis method less dependent 
of the slopes. 

From our initial results of 10 stars from automated W fitting, we show the 
age-metallicity relationship (AMR) in Fig. [5J Ed vardsson et all (|l993fi reported 
a real scatter of metallicity with respect to age, and much published literature 
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Figure 5. Age-metallicity relationship based on the 10 subgiants (solid 
points). Age uncertainty larger than 5 Gyr is shown as an arrow. Sun is 
shown as solar symbol for reference. 



states no tight correlation in AMR for [Fe /H] . Our present results here show no 
surprising contradictions. 



5. Future Work 

We will continue to obtain data until our targeted sample size of ~500 stars is 
achieved. We will improve our current automation algorithm for more robust 
W measurements and better accuracy in derived stellar properties. In addition 
to Fe, we will extend abundance analysis to C, O, Mg, and many other elements 
made possible by high-resolution spectroscopy. 

Once we have sufficient abundance information, we will be able to incorpo- 
rate this into a study of the star formation history of the solar neighborhood, 
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using observed values to constrain the models. We will also investigate the 
abundance distribution of heavy elements and the age-metallicity relation. 

By selecting a volume-limited sample, we are biased towards thin-disk stars 
that have smaller velocity dispersion from the plane. We are also biased towards 
more metal-rich stars due to a lower-limit cutoff in the color. In the sample 
of subgiants, we are biased towards less massive stars as they have a longer 
timescale in that evolutionary stage than their more massive counterparts. We 
will attempt to correct for these biases in our analysis. 

A useful by-product of this project will be a database of high resolution 
spectra of volume-limited sample of solar neighborhood stars, yielding a catalog 
of fundamental stellar properties, kinematics, and detailed abundances. 
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